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Abstract 

In an attempt to clarify the nature of the crossover from a Peierls band insulator to a Mott 
Hubbard insulator, we analyze ground-state and spectral properties of the one-dimensional half- 
filled Holstein Hubbard model using exact diagonalization techniques. 
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In a wide range of quasi-one-dimensional ma- 
terials, such as MX chains, conjugated poly- 
mers or ferroelectric perovskites, the itineracy 
of the electrons strongly competes with electron- 
electron and electron-phonon (EP) interactions, 
which tend to localize the charge carriers by es- 
tablishing spin-density-wave and charge-density- 
wave ground states, respectively. Hence, at half- 
filling, Peierls (PI) or Mott (MI) insulating phases 
are energetically favored over the metallic state. 
An interesting and still controversial question is 
whether or not only one quantum critical point 
separates the PI and MI phases at T = 0]. 
Furthermore, how is the crossover modified when 
phonon dynamical effects, which are known to be 
of particular importance in low-dimensional ma- 
terials ^, |j , are taken into account? 

The paradigm in studies of this subject is the 
half-filled Holstein-Hubbard model (HHM), de- 
fined by the Hamiltonian 

H= -i^(4 CT c J+lCT + H.c.) + [/^n iT n a 
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Here cj CT creates a spin-cr electron at Wannier 
site i {rii M — c\ a c ia ), b\ creates a local phonon 
of frequency uiq, t denotes the hopping integral, 
U is the on-site Hubbard repulsion, g is a measure 
of the EP coupling strength, and the summation 
over i extends over a periodic chain of N sites. 
Applying basically exact numerical meth- 



ods M, we are able to diagonalize the HHM on 
finite chains, preserving the full dynamics of the 
phonons. In order to characterize the ground- 
state and spectral properties of the HHM in dif- 
ferent parameter regimes, we have calculated the 
charge- and spin-structure factors at q = n 

S ^ = iEf^)'^^-!) (2) 
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(Sf = (riij — n,jJ/2), the local magnetic mo- 
ment L — 3((5,f) 2 ), the kinetic energy Eki n = 

-t(J2i,A c ia c i+ia + H.c.)), and the incoherent 
part of the optical conductivity 
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where j = -ietJ2 ia ( c L c i+i a - 4 
typical results are shown in Figs. 1 and 2 
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Our conclusions can be summarized as follows: 

(i) At U = the ground-state is a Peierls 
distorted state in the adiabatic limit luq — > 
for any finite EP coupling. As in the Holstein 
model of spinless fermions || , at ujo > quantum 
phonon fluctuations destroy the Peierls instabil- 
ity for small EP interaction strength g j| . Above 
a critical threshold g c (ujo), the HHM describes 
a PI with gapped spin and charge excitations. 
In the non-adiabatic strong EP coupling regime, 
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Figure 1: Staggered charge- and spin-density cor- 
relations (upper panel), kinetic energy and local 
magnetic moment (lower panel) in the ground 
state of the Holstein-Hubbard model (u>o /t = 1; 
N = 8). Results are shown at g 2 uJo/t — (dia- 
monds), 0.5 (circles), and 2.0 (squares). 



the system is typified by a charge-ordered bipo- 
laronic insulator rather than a traditional Peierls 
band insulator. The PI regime is characterized 
by a large (small) charge (spin) structure factor, 
a strongly reduced kinetic energy, and an opti- 
cal response that is dominated by multiphonon 
absorption and emission processes. 

(ii) Increasing U at fixed g, the Peierls dimer- 
ization and the concomitant charge order are sup- 
pressed. Accordingly the system evolves from the 
PI to the MI regime. From our numerical data we 
found evidence for only one critical point U c (cf., 
e.g., the development of the optical gap in the 
conductivity spectra shown Fig. 2). At U c , in our 
finite system a site parity change of the ground 
state takes place from P = +1 (PI) to P = — 1 
(MI), and both the spin and the charge gaps are 
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Figure 2: Optical absorption in the HHM (u>o/t=l; 
N = 8). Dashed lines give the integrated spectral 
weights S re9 (uj) = fg a TeB (u')du}' [normalized by 
S re9 (oo)]. 
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Table 1: Parity of the ground state of the HHM, 
Pfyo) = ^fyo), where the site inversion symme- 
try operator P is defined by Pc-^P^ 
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for i = 0, 1, . . . , N — 1. Charge and spin gaps, 
defined by A C (A^) = P (f + 1. f ) + E o( N 
1, f ) - 2P (f , f ) and A S (N) = P (f 
1) - E (f , f ), respectively, where E (N r , N±) 
denotes the ground-state energy of the system 
with TVj spin- up and N± spin-down electrons. 
Note that A c incorporates ground-state lattice 
relaxation effects and therefore differs from the 
optical gap In the infinite system, we ex- 
pect that A c = A s > for U < U c , whereas 
A c > A s = for U > U c . Results are given at 
g 2 uj = 2, uj /t = 1. 



expected to vanish in the thermodynamic limit 
(cf. Tab. I). 

(iii) Above U c , in the MI phase, the low- 
energy physics of the system is governed by gap- 
less spin and massive charge excitations. In the 



Mott Hubbard insulating regime the optical gap 
is by its nature a correlation gap. It is rapidly 
destroyed by doping the system away from half 
filling §. 
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